Rapid Measurement of Quantum Systems using Feedback Control 
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We introduce a feedback control algorithm that increases the speed at which a measurement 
extracts information about a d-dimensional system by a factor that scales as . Generalizing 
this algorithm, we apply it to a register of n qubits and show an improvement 0{n). We derive 
analytical bounds on the benefit provided by the feedback and perform simulations that confirm 
that this speedup is achieved. 
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Recently it has been shown that it is possible to in- 
crease the speed at which a measurement purifies the 
state of a quantum system by using real-time feedback 
control as the measurement proceeds [1 121 El HI |5l [6[ [3 
[HI [9]. Specifically, by using feedback to keep the state 
of the system diagonal in a basis that is unbiased with 
respect to the measured observable, one can make the 
system purity increase deterministically at a rate faster 
than the increase of the average purity by measurement 
alone [U [2]. Since this protocol requires an unbiased 
basis, it exploits a purely quantum mechanical effect. 
However, as a consequence this quantum feedback pro- 
tocol prevents the observer from obtaining full informa- 
tion about the initial preparation of the system, and is 
therefore not appropriate for use in applications such as 
communication channels [TO]. It is for this reason that 
the effect of the protocol is termed rapid-purification and 
not rapid-measurement 

Here we present a new protocol that can be applied to 
both quantum and classical systems, and that, in con- 
trast to quantum rapid-purification, increases the rate 
at which the observer gains information about the ini- 
tial preparation, as well as the rate at which the state 
is purified. This protocol therefore achieves not merely 
rapid purification, but also rapid measurement, and can 
be used in communication channels, state stabilisation, 
read out and error correction in quantum computers. It 
is also distinct from the rapid state-discrimination pro- 
tocol introduced recently in [iP, which is only applicable 
to non-orthogonal states and thus to quantum systems. 
Our protocol can be applied to all systems with dimen- 
sion d > 2. Asumming the measured observable scales 
with we show that our protocol increases the speed of 
a measurement by a factor 0{d'^) over an unaided mea- 
surement on a qudit. This is in contrast to the 0{d) 
speed-up acheived by previous protocols [2]. We gener- 
alise our protocol for a register of n qubits, each being 
measured independently and continuously, and obtain an 
improvement 0(n). 

The evolution of the state p of a system subject to a 
continuous measurement of an observable X is given by 



the stochastic master equation (SME) [131 El 

dp = 2-fdtV[X]p^ y/¥fdWn[X]p, (1) 

where V[A]p = ApA^ - + pA^A), and n[A]p = 

Ap-\- pA^ — Tr [{A^ + A)p\ p. The measurement strength^ 
7, determines the rate at which information is extracted, 
and thus the rate at which the system is projected onto 
a single eigenstate of X [12 . We denote the continuous 
measurement record obtained by the observer as r(t), 
and dr = ^^^{X)dt + dW. We assume that the mea- 
surement strength is much smaller than the strength of 
the Hamiltonians that can be controlled by feedback, so 
that applied unitaries can be treated as instantaneous. 

For a continuous measurement, without feedback, 
there must be no degenerate eigenvalues in the observable 
in order to guarantee purification [12j. In what follows 
we will take X to have equispaced, and hence nondegen- 
erate, eigenvalues (although we note that, with feedback, 
in order to purify an arbitrary initial state, the observ- 
able need possess only two nondegenerate eigenvalues). 
This is a fairly canonical choice, as it applies to observ- 
ables such as the components of angular momentum and 
the energy of a harmonic oscillator. This means that, 
since the transformation X ^ X -\- al for real a leaves 
the SME invariant, we can always take the d eigenval- 
ues {x} of X, to be {0, 1, . . . , (i — 1}, with corresponding 
eigenstates \x). We will label the eigenvalues of p as 
in decreasing order: Aq > Ai > • • • A^^-i. We will assume 
that p is initially completely mixed, and that any feed- 
back merely acts so as to swap states in the eigenbasis 
{\x)} of X. In this case p will remain diagonal in this 
eigenbasis, and we will denote the eigenstate of p with 
eigenvalue Xj as \xj). 

The rapid-measurement algorithms we present are 
based on the following two observations: 1. The rate at 
which a measurement distinguishes between two eigen- 
states, \x) and \x') is proportional to ^{x — x')^ . 2. As 
the measurement proceeds, Aq approaches unity. With 
no feedback, the above facts mean that at long times the 
next-largest eigenvalue Ai will be associated with adja- 
cent eigenvalues of X\ xi = ± 1 [2 . However, fact 1 



2 



means that if we apply a unitary operation to the state 
matrix so as to rearrange the eigenbasis, we can make 
\xi — xo\ as large as d — 1, leading to a — 1)^ increase 
in the rate at which we distinguish between the most and 
next-most likely state. This is the basic idea behind our 
rapid measurement algorithms. 

Previous analyses of rapid purification have considered 
both the impurity L = 1 — [1, 2 and the infidelity 

A = 1 — Ao [3 . In the long-time limit these are propor- 
tional {L ~ 2 A <^ 1) and for the most part we use A 
rather than L in this Letter. As well as considering the 
time r required for the average infidelity (A(r)) to reach 
a specified level e, one can also consider the average time 
(T) for a system to attain A(T) = e [3 . In many circum- 
stances the latter is the more useful quantity [3]. It is also 
much easier to calculate numerically. The reason is that 
the time T at which A(T) = e has well-behaved statistics 
[3], in contrast with A(t) which has extremely long tails 
at relatively large values of A. That is, the average (A) 
is greatly infiuenced by the rare cases that are slow to 
purify. Because of this, there is substantial disagreement 
between r and (T) for a qubit. It was shown in Ref. [3 , 
however, that good agreement is found between (T) and 
r', defined as the time required for (InA(r')) to reach a 
specified level In e. This is because taking the logarithm 
de-emphasizes the tails, and indeed for a qubit In A has 
a near-normal distribution [3 . 

In this paper most of our calculations concern the av- 
erage time (T) to obtain a measurement fidelity A(T) = 
e <C 1. In addition to the reasons given above, this is 
because we can obtain the scaling of (T) from analyti- 
cal bounds for the long-time behaviour of (In A). While 
we are not able to obtain bounds on (A) or (L), we in- 
clude numerical simulations of {L) for small systems. As 
we will see, these simulations suggest that all measures 
show similar increases in measurement rate under our 
algorithm. 

We first obtain analytical expressions for (In A) for 
the case of no feedback (nfb). In the limit of interest, 
where one eigenvalue is very close to one, we can use the 
fact that In A ~ Tr[ln(l — p)]. The latter expression has 
the advantage that one can obtain an exact integral for 
(Tr[ln(l — /))]), using linear quantum trajectory theory 
[lll[T5] as in Ref. [2 . For t 7"^ this integral gives 

(lnA)nfb--47t. (2) 

From this relation and the argument given above (dis- 
cussed in more detail in Ref. [3]) we expect that the mean 
time to attain A = e is, for ln(e~^) ^ 1, 

(T)„fb = (l/47)ln(e-i). (3) 

We now turn to our feedback algorithm. As explained 
above, the basic idea is to make the eigenvalues of the 
measured quantity X as different as possible for the 
eigenstates of p with the largest and second-largest eigen- 
values. Thus we take the action of the feedback unitaries 



to be such as to make xq the smallest and xi the largest 
eigenvalue of X. That is, xq = and Xi = d — 1. Then 
from the SME ^ we calculate 

d(lnA) = -47df^^^^^f^, (4) 

where (X) = Tr[Xp]. For a given A, the measure- 
ment rate is maximized by using the feedback to max- 
imise {Xy. This is achieved by using feedback to or- 
der the eigenvalues of p as Aq, A^^-i, A^i_2, • • • ,A2,Ai. 
Here the corresponding eigenstates are X-eigenstates 
|0),|1),|2)--. ,|^i-2),|d-l). We win can this L- 
ordering [17]. The feedback algorithm that maintains 
this ordering is Locally Optimal in time (LO) in the sense 
that at any point in time L-ordering will maximise the 
expected decrease in the log-infidelity subsequent to a 
weak measurent of infinitesimal duration dt. It is worth 
noting that the feedback, applying a conditional unitary 
to the system, is equivalent to changing the measured 
observable at each time step. 

We can derive bounds on the performance of LO feed- 
back in the long-time limit from bounds on 

d-l 

(X) = ^A,(d-j), (5) 

with the constraints Yl^Zl = ^ Xj > Xk > for 
j < k. It is easy to prove that A{d/2) < (X) < A(c^- 1). 
Thus in the long-time limit where A <C 1 we have 

d{\nA)i^o = -4-fdtSi^o, (6) 
{d/2f <Si.o<{d-lf. (7) 

Using this to infer approximate bounds on (T)lo7 as ex- 
plained above, we see from Eq. (|6| that ^lo is precisely 
the speed-up-factor that this feedback algorithm gives 
over the measurement without feedback: 

{Tko = {T)nfh/SLO- (8) 

The above argument thus predicts that the achievable 
speed-up in (T) increases quadratically in the size of the 
system. This scaling, which we confirm with numerical 
simulations below, is the main result of this Letter. 

We expect the true speed-up factor S'lo to be closer 
to the lower bound in Eq. Q, for the following reason. 
For A <C 1, the observer is almost certain that the sys- 
tem is in state |0), but it could be in one of the states 
• • • , |(i — 1). For some state \x) with x > 0, the larger 
X is, the better the measurement is at determining that it 
is not in that state. Hence the LO feedback ensures that 
the measurement will tend to reduce Xj-i more than Aj, 
Vj > 0. Since Aj_i > Aj, the measurement plus LO feed- 
back will tend to equalize the Xj for j > 0. The lower 
bound is attained when they are all equal. 
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FIG. 1: (color online). The speed-up in the mean time to 
reach a given infidelity for systems of dimension d — 3,4 and 
5, when the initial state is completely mixed. The dotted lines 
give the numerically calculated asymptotic speedup. 
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FIG. 2: (color online). The asymptotic speed-up for reaching 
a given level of infedelity as a function of the system dimen- 
sion. The dotted line is a linear fit to the data points. The 
solid lines are the upper and lower bounds derived in the text. 

We now compare the above calculation to numerics. 
We calculated (T)^^ and {T)j^q by stochastic simula- 
tions. The former was found to agree with the theory 
above Eq. ([3|. In both cases, we calculated the small-e 
asymptotic speed-up by extrapolating the numerical val- 
ues of (T). In Fig. [l] we plot the speed-up factor from 
Eq. ([8| as a function of e, for d = 3,4,5, also with the 
asymptotic values. This shows that most of the speed- 
up is achieved for e = 10~^, although the gap appears to 
grow with d. In Fig.[2]we plot the asymptotic speed-up as 
a function of d. As expected, it fits within the bounds of 
Eq. ([7|, and confirms a very nearly quadratic speed-up: 
the fit shown is 6' = 0.45^^ - OM + 1. 

We also consider an alternative definition for the speed- 
up, based on the time taken to achieve a fixed average 
impurity (L). We calculate this speed-up factor for the 
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FIG. 3: (color online). The asymptotic speed-up in the time 
required to reach a given level of average purity, as a function 
of the final average linear entropy. Solid lines: continual feed- 
back for systems of dimension 3, 4 and 5; Dotted line: periodic 
feedback at intervals of t = O.2/7. 

LO algorithm numerically, and plot it as a function of 
(L) in Fig. |3l While our run time is necessarily limited, 
these results indicate that the asymptotic speedup also 
falls within the bounds derived for (In A) above. 

In an experimental implementation it may not be fea- 
sible to continually, and thus very rapidly, rearrange the 
elements of p. It turns out that rearranging the ele- 
ments of p at regular intervals instead is sufficient to 
obtain most of the speed-up. In this alternative strat- 
egy the uncertainty tends to be localized into three adja- 
cent eigenstates by allowing the measurement to proceed 
without feedback for a fixed time. After the fixed time 
the feedback algorithm re-arranges the eigenvalues into 
L-ordering. We simulated this scenario for d = 4, with 
re-arrangements performed at intervals oft = O.2/7. The 
resulting speed-up is plotted in Fig. [3) showing that the 
performance is not reduced significantly by this change. 
We find that when the feedback is turned on the rate of 
decrease of (L) is in fact greater than that for the LO 
feedback algorithm. However, this rate decreases over 
time, and at no time does the alternative strategy beat 
the locally optimal algorithm. This suggests that the 
latter may be globally optimal. 

We now generalize the above rapid measurement pro- 
tocol to a register of n qubits, where each qubit is in- 
dependently and continuously measured. Instead of one 
observable X, we now have n, given by Zr = 0/^^^ ^ 
. . . ai"^^ . . . (g)/(^), where r labels the rth qubit. The SME 
describing such a measurment is [18] 

dp = ^2ndtV [Zr] p + dWrH [Zr] p. (9) 

r 

We transform Zr according io Zr ^ Zr — I and order the 
eigenvalues of p such that Aq, > A/3 when a < where 
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FIG. 4: (color online). The asymptotic speed-up in the mean 
time for a quantum register to reach a given level of infidelity, 
as a function of the number of qubits in the register. The 
dashed line is a linear fit. The solid lines are the bounds de- 
rived in the text. Inset: the optimal eigenvalue arrangement 
for a) a two qubit register, and b) a three qubit register. 



a, /3 are binary strings. Without loss of generality we can 
assume the largest eigenvalue of p (Aq) to be placed such 
that it corresponds to |0) = |00...0). Using a similar 
technique to that described above Eq. ([s]), we obtain a 
long-time analytical expression for (In A) in the absence 
of feedback: (ln(A))nfb ^ —Wtzt. For the LO feedback 
algorithm, the SME ([9| is used to calculate the average 
rate of change of the log- infidelity, which is 



time limit (A << 1) we thus find 

d{\nA)i^o = -4i^dtSi^o, (11) 
712^^/4(2^ - 1)2 < < n. (12) 

For n >7 the lower bound is approximately n/4. These 
are rigorous bounds on the speed-up factor for (In A), 
and once again we expect these to well-approximate the 
behavior of (T) . To confirm this we performed numerical 
simulations of the nfb and LO algorithms for quantum 
registers of different sizes, and extrapolated the asymp- 
totic speed-ups as above. The results are displayed in 
Fig. [4j and show the expected linear dependence on the 
number of qubits; the linear fit is = 0.718n. 

In summary, we have shown that it is possible to use 
feedback control during a measurement process to in- 
crease the speed of the measurement by a factor pro- 
portional to (P. One can also use this method to speed 
up the measurement of a register of n qubits by a factor 
that scales as n. We expect that approximate algorithms 
which are experimentally feasible should achieve most 
of this speedup, so that feedback-enhanced measurement 
could be applied in a variety of quantum technologies. 
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Now we wish to maximize the average reduction of 
the log- infidelity for a given p using feedback. This 
is acheived by reordering the elements of p, so as to 
maximize X]^(^r)^, in the following way. By defini- 
tion (above) the largest eigenvalue is at |0). The second 
largest eigenvalue A00...01 is then placed at |I) such that 
it is the maximum Hamming distance ^6j away. The 
next n largest eigenvalues are placed at one Hamming 
unit away from |1). Then the next ^€2 largest eigen- 
values are placed two Hamming units away from |1) and 
so on (this process is repeated a total of n — 1 times). 
We call this ordering H-ordering [19 and it is LO for 
the (In A) measure. Example H-orderings for a two- and 
three-qubit register are depicted in the Fig. [4] insets. 

We now bound, for a register of qubits, the amount 
by which the H-ordering algorithm speeds up the mea- 
surement process. To do this we must bound X]^(^r)^- 
The upper bound corresponds to the probablity being 
concentrated into the eigenvalue that is maximally dis- 
tant from |0), that is A00...1 = A, while the lower bound 
is reached when the probablity A is distributed equally 
over the remaining 2"^ — 1 eigenvalues. The bounds are 
(n22^/(2^ - 1)2)A2 < Er(^r)^ < 4nA2. In the long- 
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